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CORRELATION OF ELASTOHYDRODYNAMIC FILM THICKNESS 
MEASUREMENTS FOR FLUOROCARBON, TYPE II ESTER, 

AND POLYPHENYL ETHER LUBRICANTS 
by Stuart H. Loewenthal, Richard J. Parker, and Erwin V. Za ret sky 

Lewis Research Center and 
U. S. Army Air Mobility R&D Laboratory 

SUMMARY 

A minimum film thickness correlation for heavily loaded elastohydrodynamic (EHD) 
contacts was formulated from X-ray film thickness measurements made with fluoro- 
carbon, type II ester, polyphenyl ether, and synthetic paraffinic lubricants. The semi- 
empirical model was developed from test data which covered a wide and practical range 

Q 

of operating conditions. Maximum Hertz stresses ranged from 1.04x10 to 
2. 42x10® newtons per square meter (150 000 to 350 000 psi), disk temperatures from 
339 to 589 K (150° to 600° F), and mean surface rolling speeds from 9. 4 to 37. 6 meters 
per second (370 to 1480 in. /sec). 

A comparison of the predicted values of the minimum film thickness from the de- 
duced EHD formula with the X-ray results revealed that the film thickness correlation, 
despite its simplified form, represented the measured data reasonably well through- 
out the full range of test conditions. 

At elevated temperatures the fluorocarbon appeared to be the best film former fol- 
lowed by the synthetic paraffinic, polyphenyl ether, and type II ester fluids in diminish- 
ing order of film forming capability. 

Comparisons were made with a representative isothermal EHD analysis with and 
without thermal corrections for shear heating at the inlet of the contact zone. This 
comparison reconfirmed that theoretical consideration of heating effects would not ap- 
preciably influence the weak relation of isothermal film thickness to contact stress nor 
significantly diminish the magnitude discrepancy between theory and the X-ray data ex- 
cept at the more severe heating conditions. 

The semi-empirical film thickness equation contains a high- contact stress factor 
cp to correct classical film thickness formulas for the higher sensitivity of minimum 
film thickness to contact stress exhibited by the X-ray test data. This factor indicates 
that this observed higher than theoretically anticipated sensitivity is responsible for the 
predicted film thicknesses being up to twice as large as those actually measured. 

It is anticipated that the semi-empirical film thickness formula can be used with 
reasonable certainty for many rolling -element bearing and gear systems. These sys- 
tems include principally those composed of steel, employing the lubricants studied here- 
in, and those whose contact geometry approximates that upon which the model is based. 



INTRODUCTION 


The importance of maintaining a sufficient e lastohydr odynamic (E HD) film thickness 
between dynamically contacting machine elements has in recent years been more fully 
appreciated. The prediction of EHD film thickness has been the focal point of many 
theoretical and experimental investigations; it has been summarized well in refer- 
ences 1 and 2. 

The ratio of EHD minimum film thickness to composite surface roughness of the 
mating contact surfaces has become an acceptable indicator of the effectiveness of the 
lubricant film within the rolling-element contact zone. It has been shown experimentally 
that this ratio influences the fatigue life of rolling- element bearings (refs. 3 and 4). 
Predetermination of this lubricant parameter with an accurate prediction of minimum 
film thickness is of value to the designer in obtaining more realistic estimates of 
rolling-element fatigue life (ref. 5). 

The bulk of the experimental work conducted in elastohydrodynamic lubrication has 
been confined to conditions of moderate speeds up to 25. 4 meters per second (1000 in. / 

g 

sec) and moderate loads which resulted in maximum Hertz stresses to 1.24x10 newtons 
per square meter (N/m ) (180 000 psi) (refs. 6 to 9). The research of references 10 
and 11 has extended the EHD film thickness measurements to maximum Hertz stresses 
of 2. 42x10 N/m (350 000 psi) which includes the design operating range of most 
machine components such as bearings and gears. These data were obtained on a 
rolling-disk machine using an X-ray transmission technique to measure minimum film 
thickness. The film thickness measurements showed good qualitative agreement with 
full scale bearing test results (ref. 12); that is, very low film thicknesses were meas- 
ured at conditions similar to those where the bearings suffered surface damage. 

Results obtained by previous investigators showed reasonably good correlation at 
moderate speeds and loads between elastohydrodynamic theory and film thickness meas- 
urement (refs. 6, 8, and 9). In contrast, however, the data of references 10 and 11 
showed a marked deviation between predicted and experimental values of film thickness. 

In particular, at high contact stresses (i.e. , maximum Hertz stresses greater than 
Q 2 

1.38x10 N/m (200 000 psi)), the sensitivity of the film thickness to load as determined 
experimentally is far greater than that predicted by classical EHD theory of refer- 
ences 13 and 14. 

Several attempts have been made to resolve the apparent discrepancy between 
theory and experiment. A critical examination of the X-ray technique itself was made 
for possible load dependent experimental errors (ref. 15). However, no experimental 
factors were uncovered which could seriously alter the accuracy of the X-ray measure- 
ments. On the theoretical side, the influence of several possible rheological factors 
has been investigated. These factors were the effects of a non- Newtonian lubricant of 
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the Ree-Eyring form (ref. 16), the effects of heating at the inlet of the contact region 
(ref. 17), and the effects of a reduced lubricant viscosity- pressure dependence using 
both a composite exponential model (ref. 18) and a power -law model (ref. 19). 

Each of the previous modifications to elastohydrodynamic theory has succeeded 
somewhat in improving the agreement between theory and experimental data within the 
heavy load regime. However, the resulting predicted values of film thickness differed 
little in magnitude from those computed using classical EHD theory. Furthermore, the 
modified theories do not sufficiently account for the observed high film thicknesses - 
load dependence to allow accurate predictions of film thickness under realistic operating 
conditions. 

The experimental data of references 20 and 21 also show a film thickness sensitivity 

to stress greater than theoretical for maximum Hertz stresses greater than about 
Q 2 

1. 04x10 N/m (150 000 psi). These data, obtained by an optical interferometry tech- 
nique with sliding point contacts, tend to support the measurements obtained by the 
X-ray technique of references 10 and 11. 

It was the objective of the work reported herein to generalize the empirical film 
thickness model developed in reference 22 through an analysis of the experimental data 
of references 10 and 11 and to compare this derived relation with that of classical elas- 
tohydrodynamic theory. 


SYMBOLS 


a 

B 

b 


Ej,E2 


E' 


G 

min 

H* 

h . 

min 


minor semi-axis of Hertzian contact, m (in. ) 
constant, eq. (6) 

major semi-axis of Hertzian contact, m (in. ) 
coefficient, eq. (1) 

2 

modulus of elasticity of elements 1 and 2, N/m (psi) 

A - 1 - vl\ p 

1 + 2 ^ N/m^ (psi) 

\ttEi nE^J 

film thickness - stress function, eq. (2) 

n 

dimensionless speed factor, (u/N) 
nondimensional minimum film thickness, 
lubricant film thickness factor, 
minimum film thickness, m (in. ) 
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i experimental maximum Hertz stress subscript 

j test lubricant subscript 

Kj, K? ,kj lubricant coefficients in empirical film thickness formula 
N constant, 9. 4 m/sec (370 in. /sec) 

nondimensional speed - viscosity parameter exponent, eq. (1) 

J 

1 . mean value of n. ., eq. (2) 

J J 

nondimensional stress parameter, Pjjz;/®' 

2 

maximum Hertz stress, N/m (psi) 
radius of elements 1 and 2 in rolling direction, m (in. ) 




"Hz 

Phz 

Rj, R2 
R' 


U 

u 

UpUg 

a 

^ 1,2 


equivalent radius, ( ^ 


\-l 


yRl Rg^ 

disk temperature, K (°F) 


, m (in. ) 


nondimensional speed - viscosity parameter, p^u/E’R' 

mean surface velocity, l/2(uj^ + Ug), m/sec (in. /sec) 

surface velocities of elements 1 and 2, m/sec (in. /sec) 

2 - 1 

pressure- viscosity coefficient, m /N (psi“ ) 

2 2 

inlet absolute viscosity, N-sec/m , (Ib-sec/in. ) 
Poisson's ratio of elements 1 and 2 
high contact stress factor, eqs. (6) and (8) 
film thickness - stress function, eq. (4) 


X-RAY FILM THICKNESS TEST DATA 

The film thickness model presented herein was developed from the film thickness 
data (refs. 10 and 11) obtained in an X-ray rolling-disk machine. The four lubricants 
studied were a type II ester, a fluorocarbon, a polyphenyl ether, and a synthetic paraf- 
finic oil. These lubricants are considered to be promising candidates for high temper- 
ature bearing application (ref. 12). Their properties are summarized in table I. 

Briefly, the method of measuring film thickness utilizing the X-ray transmission 
technique involves flooding the contact region between the two test disks with X-rays. 
The rate of X-ray transmission which penetrates the contact within a narrow bandwidth 
is determined by a scintillation counter. Since the greatest constriction occurs at the 
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trailing edge of the contact, the X-ray count becomes a measure of the lubricant's mini- 
mum film thickness. The side lobes of the contact are avoided by the beam and thus do 
not influence the reading. 

The range of test conditions reported in references 10 and 11 include disk tempera- 
tures from 339 to 589 K (150° to 600° F), surface rolling speeds from 9. 4 to 
37. 6 meters per second (370 to 1480 in. /sec), which correspond to disk rotational 

q 

speeds from 5000 to 20 000 rpm, and maximum Hertz stresses from 1. 04x10 to 
2. 42x10^ N/m^ (150 000 to 350 000 psi). Two AISI M-50 steel disks each with a rolling 
radius of 1.83 centimeters (0. 72 in. ) and a surface finish of 2. 5x10® to 5. 0x10® centi- 
meters (1 to 2 pin. ) rms were used as the test specimens. 

Both crowned disks (with a crown radius of 27. 9 cm (11 in. )) and crowned-cone 
disks (with a cone angle of 10°, which introduce a small spin component in the contact) 
were tested, and no significant differences were reported between the two sets of film 
thickness data (refs. 10 and 11). All the data reported herein with the exception of the 
test data for the polyphenyl ether lubricant were generated with the crowned-cone test 
disks appearing in figure 1. 

FORMULATION OF EHD FILM THICKNESS CORRELATION 

The effect of surface speed u and lubricant absolute viscosity /i on the measured 
value of minimum film thickness h^^^^ at each experimental contact stress level Pjj^^ 
for the four test lubricants is shown in figure 2. The summary plot for the synthetic 
paraffinic lubricant in figure 2(d) has been presented before in reference 22. In figure 2 
the log of thfe dimensionless film thickness parameter - ^min^^' plotted 

against the log of the dimensionless speed- viscosity parameter U - (p^u/E'R'). Be- 
cause of the linear relation exhibited by the test data on this logarithmic plot as deter- 
mined by a linear regression analysis, there exists the simple power relation 


^^min^i, ; 


C. 

1,3 


( 1 ) 


where the subscript i = (1 — 5) designates one of the five experimental maximum Hertz 
stress levels, the subscript j = (1 - 4) designates one of the four test fluids, and U 
ranges from the lowest to highest experimental value. 

When comparing the plots for each lubricant in figure 2 it is apparent that both the 
polyphenyl ether and synthetic paraffinic film thickness data, unlike that for the fluoro- 
carbon or type II ester fluids, display a somewhat increased sensitivity to rolling sur- 
face speed and lubricant viscosity with increasing contact stress. This is reflected by 
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the increase in the slope of lines with higher contact stress appearing in figures 2(c) 
and (d) or more explicitly by the variation of exponent n- ■ in equation (1) which is 

J 

tabulated along with coefficient C. ■ in table II. 

— 

This variation of the dependence of to U with increasing is not anti- 

cipated from theoretical considerations. In view of the overall consistency of the test 
data, it is not apparent why this anomaly was observed for just two of the four test 
lubricants. However, the effect of neglecting such an exponential variation in the case 
of the synthetic paraffinic and polyphenyl ether fluids will have a minimum effect on the 
resulting generalized film thickness expression as will be shown later. Accordingly, 
selecting mean values of exponents n. ^ and n, . results in the values of exponent h. 
listed in table III. These values have been determined for each of the test lubricants 
from the X-ray test data. 

Equation (1) can now be written in the following form: 


(«min)j - 


U ^f(p 


Hz^j 


( 2 ) 


where the presently unknown continuous function been introduced to describe 

the dependence of minimum film thickness on the peak contact pressure for each lubri- 
cant. The main objective in this approach is to separate the effects of maximum Hertz 
stress (contact pressure) on the value of film thickness from those effects contributed 
by surface speed and lubricant viscosity. Having thus isolated the effects of contact 
pressure, there remains the task of representing the influence of contact pressure on 
film thickness by a single mathematical expression for all test fluids, that is, formulat- 
ing f(Pjj 2 )- 

By combining equations (1) and (2) an expression can be written for f(pjj^)j where 


f(p 


Hz^i 


- C, 




i 




or 



( 3 ) 


It is now possible to evaluate f(Pjjz)j the experimentally deduced values of the ex- 

pression appearing on the right side of equation (3). 

It is evident from equation (3) and the previous discussion that f(pjj^)j is an ex- 
plicit function of U. The extent of this dependence is illustrated in figure 3 which shows 
the maximum variation of f(Pjj^) with evaluated at both the maximum and mini- 

mum experimental value of U for each of the test fluids. In the case of the fluorocar- 
bon and type II ester fluids, the effect of U on the f(pjj^)j term in equation (2) is neg- 
ligible as evidenced by the nearly coincident lines in figures 3(a) and (b). With regard 
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to the polyphenyl ether and synthetic paraffinic test fluids, the sensitivity of minimum 
film thickness to maximum Hertz stress does vary somewhat with changes in operating 
surface speed and lubricant viscosity as evidenced by figures 3(c) and (d). This varia- 
tion is anticipated since it stems from the variation of the U exponent n^ ^ with Pjj^j 
which was discussed earlier. The influence of this variation on the accuracy of the re- 
sulting film thickness expression does not warrant the introduction of an additional U 
dependent factor into the finalized correlation as will be shown later. 

Quite satisfactory results can be obtained by evaluating the f(pjj^) in equation (3) at 
the mean experimental values of U. This has been accomplished in figure 4 where the 
relation of f(Pjj 2 )j) which has been evaluated at Uj^^g^n’ ^Hz is shown. It is ap- 
parent that there is a great similarity in the shape of the curves appearing in this figure . 
That is to say, the effect of maximum Hertz stress on the ratio of (H • ). to 

U •’ is nearly the same, apart from some constant multiplier, say k., for all the fluids 
tested. Thus, the effect of contact pressure on the test lubricants minimum film thick- 
ness can be represented by a single generalized function xjy by the relation 

o 




^^Phz^] 


(4) 


where constant k. is an arbitrary lubricant constant used to normalize the value of 
f(PH 2 )j at =3. 09X10"^ - that is, at Pjj^ = 1.04x10^ N/m^ (150 000 psi). Table III 
Lists the respective kj constants for each lubricant. 

Having a function i//g which satisfactorily describes the effects of Pj^^^ on 
the generalized film thickness expression can be completed by substituting equation (4) 
into equation (2): 




(5) 


This expression is, in itself, a film thickness correlation which can be used to satis- 
factorily forecast minimum film thickness at high contact stress levels. However, as 
a matter of user convenience, equation (5) will be altered slightly to a more universal 
form. 


FORMULATION OF EHD HIGH- CONTACT- S TRESS FACTOR 

As previously discussed, film thickness values forecasted by currently accepted 
EHD theory have shown reasonably good agreement with experimental data for maximum 
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Q 2 

Hertz pressures less than approximately 1. 04x10 N/m (150 000 psi). Above this 
stress level the conventional theory seriously overestimates the extent of the film gen- 
erated by the lubricant as evidenced by test data (refs. 10 and 11). Thus, it is most 
desirable to introduce some factor to adjust current film thickness formulas for the 
deviation between theory and experiment at high applied loads. 

It is generally recognized that film thickness is only moderately dependent on con- 
tact stress at the lower stress levels. Typically, for line contact the film thickness is 
proportional to maximum Hertz stress to the -0. 22 power , that is, H oc > 

where the stress parameter exponent arbitrarily selected here comes from the iso- 
thermal theory of Cheng (ref« 23). This proportionality can be introduced into the film 

thickness relation shown in equation (5) by simply defining a new factor such that 

s 


= B 


4'. 


CP W®* 22 


( 6 ) 


o _ A 22 

where constant B has been arbitrarily chosen to equal (3. 09x10" ) to make 
— -3 

<P g = = 1 at = 3. 09x10 . Incorporating equation (6) into equation (5) yields 


min j Hz 


where K? = kj/B. Coefficient K?" together with exponent hj and factor kj are listed 

in table III. The parameter is referred to as an E HD high- contact- stress factor. 

s 

It is, in essence, a measure of the additional reduction of lubricant minimum film thick- 
ness not forecasted by commonly used film thickness formulas. 

Figure 5 shows the variation of factor (p „ with contact stress together with the 

o 

following polynomial expression which closely fits this curve; 

^s " 27. 5xlO^Pjj^) + 0.806 (8) 

It can be seen that at the lower contact stresses where film thickness sensitivity to 
contact stress has been shown to be in accordance with theory (refs. 6, 8, 9, and 24) 
the value of the cp factor is close to unity. With increasing contact stress the dis- 

D 

parity between theory and measurement becomes more significant as evidenced by the 
diminishing value of cp^. In fact, at the maximum experimental contact pressure of 
2. 4x10® N/m^ (3. 5x10^ psi), it is apparent from figure 5 that the measured value of 
film thickness is approximately half that predicted by film thickness formulas developed 
from classical theory. 
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In viewing the variation of the U exponent n., appearing in table III, it is apparent 
that the value of fij for the first three test fluids are nearly equal, averaging approxi- 
mately 0. 62. However, the value of hj for the polyphenyl ether fluid set at 0. 83 is 
significantly higher than the rest. This result qualitatively corresponds to the optical 
film thickness experiments conducted by Westlake and Cameron (ref. 25). In refer- 
ence 25 the speed- viscosity parameter exponent of a similar polyphenyl ether fluid, at 
a value of 0. 82, was found to be somewhat larger than the exponent of any other fluid 
tested including that for a fluorocarbon and synthetic paraffinic lubricant. However, it 
should be pointed out that in contrast to the results of the present work no variation of 
hj with Pjj^; foi" the minimum film thickness case was observed in reference 25 for 
either the polyphenyl ether or synthetic paraffinic oils. Because of the apparent differ- 
ences in operational conditions in terms of inlet shear rates, contact geometry, and, 
most importantly, contact stress levels, a definitive comparison between the results of 
optical film thickness measurements reported in reference 25 and those arrived at by 
using the X-ray technique is not possible. 

It was determined from numerical comparisons between the X-ray test data and 
predictions from equation (7) that, due to the commonality of the values of exponent n^ 
among the test fluids, only a small loss of accuracy would result by setting iij at a 
nominal value of 0. 62 for all four test fluids. Taking advantage of this last simplifica- 
tion, equation (7) can be written in the following final form; 


min 


KjU 


0.62p-0.22 

^Hz 




(9) 


where the lubrication parameter K? has been adjusted to to reflect the change in 
the exponent of U. Table IH lists the appropriate value of Kj for use in equation (9). 


DISCUSSION OF RESULTS 

Comparison with Test Data in View of Commonly Used Film Thickness Formulas 

The results of the present analysis are compared with the X-ray test data at several 
temperatures in figure 6. In this figure, nondimensional minimum film thickness is 
plotted as a function of maximum Hertz pressure at several rolling speeds. The numer- 
ical results of this comparison are summarized in table IV. It is evident that the de- 
rived film thickness formula (eq. (9)), although reduced to a simplified form, does 
represent the measured data quite well for all four test fluids over the full range of test 
conditions. Predicted film thicknesses are generally within 0. 05 micrometer (2 pin. ) 
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or 10 percent, which ever is greater, of those measured except for those of the syn- 
thetic paraffinic oil at high and low values of U. At these conditions limits of accuracy 
are set to ±35 percent due to, in part, the simplifications made earlier. 

Comparisons made between the X-ray test data and those film thickness values 
computed from commonly used film thickness formulas have shown that these formulas 
may seriously overestimate the extent of lubricant film actually present within the con- 
tact (refs. 18 and 22). This is illustrated in figure 7 which shows a comparison from 
reference 22 of predicted film thickness from a widely used isothermal EHD formula for 
bodies in line contact and measured film thickness for a synthetic paraffinic oil. 

At the higher rolling speeds and lower disk temperatures, hence higher lubricant 
viscosity, the effects of shear heating at the inlet of the contact zone are known to cause 
appreciable film thinning (refs. 17 and 26). Thermal effects can be taken into account 
by applying the thermal reduction cpr^ provided by Cheng (ref. 18) to the film thickness 
values computed from the isothermal formula. Correcting calculated film thickness for 
thermal effects markedly reduces the overall magnitude discrepancy between predicted 
and measured film thickness at the most severe heating test condition (fig. 7(a)) but has 
a relatively mild effect at the lower disk speeds and at any of the elevated temperature 
test conditions as can be seen, for example, from figure 7(b). Furthermore, the 
thermal reduction factor cpr^ is a very weak function of contact pressure (ref. 17). 
Consequently, heating effects would not be expected either to appreciably influence the 
relation of isothermal film thickness to applied load or explain the high film thickness - 
load dependency exhibited by the X-ray test data. 


Effect of Lubricant on Film Thickness Correlation 

The empirical factors used to formulate the present correlation have been developed 
from the experience gained with four test fluids . At present no meaningful generaliza- 
tions can be made regarding the extension of the present deduced relation to systems 
employing different lubricant types or formulations without the benefit of additional ex- 
perimental information. On the other hand, application of the film thickness correlation 
to systems utilizing the lubricants under study over similar conditions can be made with 
reasonable confidence. 

It may be apparent that the pressure -viscosity coefficient a, which is customarily 
used to characterize the film forming capabilities of a lubricant apart from the effects 
of absolute viscosity, is conspicuously absent from the present formulation. In the pre- 
sent model, the role formerly played by a has been fulfilled in part by the lubricant 
coefficient Kj. That is, a lubricant's film forming capabilities can be ascertained by 
knowing its Kj and its absolute viscosity at a given operating temperature condition. 
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In the case of a rolling element bearing, this temperature can be taken as the tempera- 
ture of the rolling element surface of interest. The good qualitative correlation which 
exists between the deduced lubricant coefficients Kj appearing in table HI and the room 
temperature values of the pressure -viscosity coefficient a listed in table I determined 
from optical film thickness measurements (refs. 25 and 27) tends to support this conten- 
tion. 

An important distinction between a and is that a is temperature dependent 
where Kj as presently defined is not. A careful examination of test data revealed that 
the effects of temperature on minimum film thickness are adequately reflected by the 
variation in absolute viscosity and that the added complication of an additional tempera- 
ture dependent variable could thus be avoided. Furthermore, the availability of perti- 
nent pressure- viscosity data at elevated temperature under the appropriate shear rate 
and pressure conditions for film thickness calculation purposes have been generally 
limited. There has been, however, increasingly more attention directed at obtaining 
these pressure- viscosity data in recent years (refs. 20, 25, and 27). In view of the 
aforementioned, the pressure- viscosity coefficient a was eliminated in the present 
film thickness model. 


Effects of Contact Geometry and Material 

Contact geometry . - The present correlation is based exclusively on measurements 
made with a single disk geometry chosen to simulate the ball inner race contact of a 
120- millimeter-bore angular- contact ball bearing (ref. 28). The contact between the 
test disks approaches the condition of line contact with an ellipticity ratio b/a of 5. 9 
where b and a are the major and minor semi-axes of the contact ellipse, respectively. 
The equivalent radius of curvature in the direction of rolling R' for the test disks is 
0. 915 centimeter (0. 36 in. ). 

It is difficult to extend the results presented herein to different contact geometries 
with complete assurance without further experimental verification. However, from a 
practical standpoint, it is anticipated that the overall effect of contact geometry on the 
value of film thickness is minimal. Cheng (ref. 23) has theoretically shown that the 
proportions of the contact ellipse with b/a varying from 1 to 5 have a relatively mild 
effect on film thickness; that is, the dependence of film thickness on the mean surface 
speed u, absolute viscosity and the contact stress Ppj^ changes little as the shape 
of the contact ellipse varies from point to line contact. Similarly, Archard and Cowking 
(ref. 29) have established that there is great similarity between the EHD lubrication of 
point and line contacts. A second factor is that the line contact geometry selected for 
the X-ray test disks does in fact simulate the contact shape that' exists between the races 
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and the balls or rollers in rolling-element bearings as well as the contact shape between 
gear teeth which usually approximate the line contact case. In view of these considera- 
tions, the semi-empirical minimum film thickness formula presented can be used with 
reasonable certainty for most practical applications without further modifications for 
small differences in contact geometry. 

With regard to the size of the contacting elements, elastohydrodynamic theory indi- 
cates that film thickness is moderately dependent on the contacting elements' equivalent 
radius of curvature in the rolling direction R* (ref. 1). The current film thickness is a 
function of R' to the 0. 38 power at a given contact stress level. This value is approx- 
imately in accordance with conventional EHD theory. However, it is recognized that the 
sensitivity of to R’ in the heavy load regime remains to be established experi- 

mentally. Until such time, minimum film thicknesses forecasted by equation (9) will be 
most successful for those systems in which the R' of the contact approximates that of 
the system under study. This limitation is not anticipated to seriously encumber the use 
of this film thickness expression for the majority of applications. 

Material. - The effects of material properties in terms of Young's modulus on film 

9/2 5 

thickness at high contact pressures (up to 3.45x10 N/m (5x10 psi) maximum Hertz 
stress) have been demonstrated by experiment (ref. 30) to be minimal. These tests 
which confirm theoretical expectations were conducted by Gohar (ref. 30) using inter- 
ferometry to measure the film generated between a rolling steel ball and a flat glass 
plate. It is anticipated that the choice of materials other than steel will not appreciably 
alter the form of equation (9), unless the elastic properties of the material of interest 
are markedly different from those of steel. 


Relative Film Forming Ability of Test Fluids 

As a design tool, the film thickness correlation developed herein can be used to 

screen prospective lubricants for particular applications once the lubricant coefficient 

K. has been experimentally ascertained. For example, the relative film forming ability 

of the test lubricants can be ranked at elevated temperatures under a practical contact 

load condition as shown in figure 8. In this plot the lubricant film thickness factor H* 

defined as the minimum film thickness H^. divided by a dimensionless speed factor 

Inin Q c\ 

G is plotted against the disk temperature at a maximum Hertz stress of 2. 07x10^ N/m^ 
(300 000 psi). It should be mentioned that the use of the lubricant film thickness factor 
H* in this figure does not alter the relative ranking of the lubricants made in this com- 
parison. This parameter serves mainly as a convenient scale factor to account for the 
beneficial effects of speed which affect all of the lubricants to the same extent. 
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It is apparent that the fluorocarbon lubricant stands out as the best film former at 
the elevated temperatures. This is attributed to the fluorocarbon fluid's high absolute 
viscosity in conjunction with the fact that it had the largest lubricant coefficient Kj of 
all the lubricants examined. 

The synthetic paraffinic oil without additives shown in this figure generates about 
one -third as much film as does the fluorocarbon lubricant. The presence of an organic 
phosphonate antiwear additive to improve boundary lubricating characteristics would 
make the synthetic paraffinic oil as good a film former as the fluorocarbon according 
to the X-ray film thickness data published in reference 10. This observation is further 
strengthened by full-scale bearing tests (ref. 12). These tests showed that a 120- 
millimeter -bore angular- contact ball bearing could achieve long term operation with 
either the synthetic paraffinic oil containing an additive or the fluorocarbon lubricant at 
temperatures to 589 K (600° F) under a comparable heavy load condition (maximum 
Hertz stress of inner race ball contact of approximately 2. 07x10 N/m (300 000 psi)). 


CONCLUDING REMARKS 

The existence of the unexpectedly high dependence of film thickness on maximum 
Hertz stress within a highly loaded concentrated contact as exhibited by the X-ray test 
data has not yet received universal acceptance. As discussed previously, the reserva- 
tions which exist concerning this phenomena stem largely from the absence of a well 
defined theoretical explanation. Apart from the experimental investigations cited pre- 
viously there is at present a scarcity of published film thickness test data under condi- 
tions of high- contact stress which would help to resolve this controversy. Nonetheless, 
the results from full-scale bearing tests conducted in reference 12 tend to affirm the 
existence of such a high load - film thickness dependence. In this investigation, a series 
of 120- millimeter-bore angular -contact ball bearings were subject to long term fatigue 
tests at elevated temperatures. These bearings had contact geometries which were 
modeled by the test disks used in the X-ray experiments (refs. 10 and 11). 

Preliminary tests were conducted with a S 3 mthetic paraffinic oil similar to that 
studied herein. The inner race had an operating temperature of 589 K (600° F) and an 
inner raceway maximum Hertz stress 2. 36x10 N/m (342 000 psi). Inspection of 
several bearings revealed that appreciable ball wear had occurred and that surface 
glazing was apparent on the races. It was concluded from this inspection that these 
bearings were operating in a mixed boundary elastohydrodynamic lubrication regime 
with boundary lubrication being the predominant mode. Similar conclusions were drawn 
when repeating this preliminary test with a fluorocarbon lubricant but at a slightly lower 
maximum Hertz stress level of 2.23x10^ N/m^ (323 000 psi). 
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As a corrective measure, the thrust load acting on these bearings was reduced to 
assure predominantly EHD lubrication for the remainder of the fatigue tests. This re- 
sulted in a reduction of the maximum Hertz stress acting on the inner race to approxi- 
mately 2. 23x10^ N/m^ (323 000 psi) and 2. 04x10^ N/m^ (295 000 psi) for the bearing 
operating with the synthetic paraffinic and fluorocarbon lubricants, respectively. 

Subsequent examination under these new load conditions of an unfailed test bearing 
run with the synthetic paraffinic oil indicated that there was no measurable wear or 
weight change of bearing components. On the basis of this post-test examination, which 
also included surface trace measurements, together with the appearance of the rolling- 
element surfaces it was concluded that almost complete EHD lubrication existed during 
the testing of these bearings. Post- test examination of the fluorocarbon lubricated 
bearings revealed that the predominant mode of lubrication had also been affected. The 
lubrication mode ranged from boundary to elastohydrodynamic with some, test bearings 
exhibiting extremely good surface appearances and others showing considerable glazing. 

The dramatic improvement in the quality of lubrication in the case of the synthetic 
paraffinic lubricated bearings and in those operated with the fluorocarbon fluid to a 
lesser extent as a result of the small reduction in the operating raceway contact stress 
level tends to give physical confirmation to the existence of a high dependency of film 
thickness to contact stress. 

In terms of classical EHD theory, the reported reductions in maximum Hertz stress 

ft 2 

levels acting on the synthetic paraffinic (1.31x10 N/m (19 000 psi)) and fluorocarbon 
8 2 

(1. 93x10 N/m (28 000 psi)) lubricated bearings would result in an increase of lubri- 
cant film of only 1.3 and 2. 0 percent, respectively. Obviously, this increase in film is 
not sufficient to account for the observed shift in the system's mode of lubrication. In 
contrast to classical theory, the contact pressure terms appearing in the heavy load 
film thickness formula (eq. (9)), that is, (Pjj 2 ) ^ gy would yield films that were 
19. 2 and 22. 0 percent greater for the synthetic paraffinic and fluorocarbon lubricants, 
respectively. These predicted increases in film thickness are undoubtedly more in line 
with those actually experienced. 


SUMMARY OF RESULTS 

A minimum film thickness correlation for heavily loaded elastohydr odynami c con- 
tacts was formulated from X-ray film thickness measurements made with synthetic par- 
affinic, fluorocarbon, type II ester, and polyphenyl ether test fluids . The film thick- 
ness test data covered a wide and practical range of operating conditions. Maximum 
Hertz stresses ranged from 1. 04x10^ to 2. 42x10^ N/m^ (150 000 to 350 000 psi), disk 
temperatures from 339 to 589 K (150^ to 600° F), and mean surface rolling speeds from 
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9. 4 to 37. 6 meters per second (370 to 1480 in. /sec). The deduced film thickness for- 
mula contained a high- contact- stress factor cp ^ to correct classical EHD formulas for 
the high sensitivity of minimum film thickness to contact stress exhibited by the test 
data under heavy loads . 

Predicted values of minimum film thickness were compared to X-ray film thickness 
measurements and contrasted against the results from a well known isothermal EHD 
analysis with and without thermal corrections for shear heating. The effects of contact 
geometry, material, and lubricant properties on predicted film thickness were consid- 
ered. The following results were obtained: 

1. The film thickness correlation represented the X-ray test data reasonably well 
throughout the full range of test conditions. 

2. The observed greater dependency of film thickness to maximum Hertz stress 
than that accounted for by classical EHD theory is responsible for film thicknesses 
computed by commonly used formulas to be up to twice larger than those actually meas- 
ured. Corrections can be made by applying the EHD high- contact stress factor cp 
developed herein to these computed values. 

3. It is anticipated that the empirical film thickness formula can be used to fore- 
cast minimum film thickness under heavy loads with reasonable certainty for many 
rolling-element bearing and gear systems. These include principally those systems 
employing the lubricants studied herein and whose contact geometry approximates that 
on which the model is based. 

4. A theoretical consideration of lubricant shear heating at the inlet of the contact 
zone does not appreciably alter the weak relation of isothermal film thickness to contact 
stress nor significantly diminish the magnitude discrepancy between theory and the 
X-ray data except at the most severe heating conditions. 

5. At elevated lubricant temperatures the fluorocarbon lubricant appears to be the 
best film former followed by the synthetic paraffinic (without an antiwear additive), 
polyphenyl ether, and type II ester fluids in diminishing film forming order. 

Lewis Research Center, 

National Aeronautics and Space Administration, 
and 

U. S. Army Air Mobility R&D Laboratory, 

Cleveland, Ohio, August 5, 1974, 

501-24. 
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TABLE I. - PROPERTIES OF TEST RUBRIC ANTS®^ 


Property 

Lubricant designation 


T 3 rpe II ester 

Fluorocarbon 

Polyphenyl ether 

Synthetic paraffinic oil 

Additives 

Oxidation inhibitor; 
corrosion inhibitor; 
antiwear additive 

None 

Oxidation inhibitor 

None 

Kinematic viscosity, 
cS (or 10'® m^/sec): 





At 233 K (-40° F) 

8900 

>100 000 

>100 000 

>100 000 

At 311 K (100° F) 

29.0 

298 

358 

443 

At 372 K (210° F) 

5.4 

29.8 

13.0 

39. 7 

Flash point, K (°F) 

533 (500) 

None 

555 (500) 

542 (515) 

Autoignition tempera- 
ture, K (°F) 

716 (830) 

None 

886 (1135) 

703 (805) 

Specific gravity at 
478 K (400° F) 

0.85 

1.59 

1. 05 

0.74 

Pressure- viscosity 
coefficient at 298 K 

^1.6x10'® (1. IxlO'^) 

°3.6xl0'® (2.5X10'^) 

®2.8xl0~® (1.9X10'^) 

’^1.6x10'® (1. IxlO'^) 

(77° F), m^N (psi'^) 






^From ref. 11. 
'^From ref. 27. 
‘^From ref. 25. 
























TABLE n. - EMPIRICAL MINIMUM FILM THICKNESS COEFFICIENTS FOR EQUATION (1) 


Nondimensional 
stress parameter, 

^Hz 

Maximum Hertz stress, 
^Hz 

Synthetic 

paraffinic 

Fluorocarbon 

Type II ester 

Polyphenyl ether 


"i2 

% 

"^3 

^<4 

"^4 

N/m^ 

psi 


•^il 

3. 09x10"^ 

1. 04x10® 

1.5X10® 

28. 1 

0.58 

90.6 

0.60 

35.8 

0. 60 

58. 1 


4. 12 

1.38 

2.0 

44. 6 

. 61 

109. 

.61 

34.2 

.60 

7. 50X10® 

.81 

5. 14 

1. 72 

2. 5 

93.5 

.65 

127. 

.62 

36.8 

.61 

5.25X10^ 

.90 

6. 17 

2. 07 

3. 0 

231. 

.70 

64.6 

.61 

28. 1 

.61 

5. 83x10® 

1.0 

6.69 

2.24 

3.25 


— 

38.8 

.59 

— 

— 


— 

7.20 

2.42 

3.5 

655. 

.76 


— 

12. 6 

.59 


— 


TABLE m. - MEAN VALUES OF THE DIMENSIONLESS SPEED VISCOSITY 
PARAMETER EXPONENT AND LUBRICANT 
PARAMETERS k., K? , AND K. 


Parameter 

Type II ester 

Fluorocarbon 

Polyphenyl ether 

Synthetic paraffinic 

5 

"j 

38. 

10. 7 
.60 
18.2 

109. 
30. 6 
. 61 
44. 8 

1. 4x10^ 
3940. 

.83 

24.9 

156. 

43.8 

.66 

18.2 


20 



























TABLE IV. - COMPARISON BETWEEN PREDICTED AND MEASURED MINIMUM FILM THICKNESS* 


(a) Synthetic paraffinic oil 


(b) Type n ester lubricant 


Disk speed, 
rpm 

Maximum Hertz stress 

Disk temperature, K (°F) 

N/m^ 

psi 

339 (150) 

422 (300) 

478 (400) 

Nondimensional minimum film thickness, H • 

min 

X-ray 

data 

Equa- 
tion (9) 

X-ray 

data 

Equa- 
tion (9) 

X-ray 

data 

Equa- 
tion (9) 

5 000 

1. 04X10® 

1.5x10® 

98X10'® 

79x10'® 

21x10'® 

19X10'® 

14x10'® 

11x10'® 


1.72 

2.5 

80 

60 

13 

14 

8 

9 


2.4 

3.5 

46 

30 

4 

7 

2 

4 

10 000 

1.04 

1.5 

U6 

122 

36 

29 

25 

17 


1.72 

2.5 

96 

92 

24 

22 

19 

13 


2.4 

3.5 

55 

47 

11 

11 

7 

7 

20 000 

1.04 

1.5 

142 

186 

48 

45 

39 

27 


1.72 

2.5 

110 

140 

35 

34 

25 

20 


2.4 

3.5 

67 

71 

17 

17 

16 

10 


Disk speed, 
rpm 

Ma.ximum Hertz stress 


Disk temperature, K 

°F) 


N/m^ 

psi 

339 (150) 

366 (200) 

422 (300) 

Nondimensional minimum film thickness, 

Hmin 

X-ray 

data 

Equa- 
tion (9) 

X-ray 

data 

Equa- 
tion (9) 

X-ray 

data 

Equa- 
tion (9) 

5 000 

1, 04X10® 

1.5x10® 

25x10'® 

26x10'® 

12x10'® 

to 

X 

o 

7x10'® 

7x10'® 


1.72 

2.5 

19 

20 

9 

9 

5 

6 


2.07 

3. 0 

U 

10 

5 

5 

3 

3 


1.04 

1.5 

42 

40 

22 

19 

12 

11 


1.72 

2.5 

31 

30 

16 

14 

9 

9 


2.07 


18 

15 

8 

7 

5 

4 



1.5 

53 

61 

29 

29 

18 

17 


1.72 

2.5 

41 

46 

25 

22 

15 

13 


2.07 

3.0 

22 

23 

15 

11 

7 

7 


(c) Fluorocarbon lubricant (d) Polyphenyl ether lubricant 


Disk speed, 
rpm 

Maximum Hertz stress 

Disk temperature, K (°F)' 

N/m^ 

psi 

422 (300) 

478 (400) 

534 (500) 

5S9 (600) 


Nondimensional minimum film thickness, H 

min 


X-ray 

data 

Equa- 
tion (9) 

X-ray 

data 

Equa- 
tion (9) 

X-ray 

data 

Equa- 
tion (9) 

X-ray 

data 

Equa- 
tion (9) 

5 000 

1. 04x10® 

1.5x10® 

68x10'® 

60X10'® 

31X10'® 

35X10'® 

21X10'® 

23X10'® 

16x10'® 

16X10'® 


1.72 

2.5 

46 

45 

20 

26 

17 

17 

12 

12 


2.07 

3.0 

31 

35 

15 

20 

13 

13 

9 

10 

10 000 

1.04 

1.5 

90 

92 

49 

53 

38 

35 

26 

25 


1.72 

2.5 

73 

70 

36 

40 

31 

27 

21 

19 


2.07 

3.0 

53 

54 

26 

31 

25 

21 

15 

15 

20 000 

1.04 

1.5 

125 

142 

79 

82 

56 

54 

38 

39 


1.72 

2.5 

101 

108 

68 

62 

46 

41 

31 

29 


2.07 

3.0 

75 

83 

48 

48 

35 

32 

25 

23 


a 


Disk speed, 
rpm 

Maximum Hertz stress 

Disk temperature, K (°F) 

N/m® 

psi 

422 (300) 

505 (450) 

Nondimensional minimum film thickness, 

x-ray 

data 

Equa- 
tion (9) 

X-ray 

data 

Equa- 
tion (9) 

5 000 

1. 04x10® 

1.5X10® 

16x10'® 

16X10'® 




1.72 

2.5 

11 

9 

4X10'® 

6X10'® 


2.4 

3.5 

8 

6 

3 

3 

10 000 

1.04 

1.5 

25 

22 




1.72 

2.5 

21 

14 

6 

10 


2.4 

3.5 

17 

9 

4 

5 

20 000 

1.04 

1.5 

36 

34 




1.72 

2.5 

32 

22 

12 

15 


2.4 

3.5 

27 

15 


7 


From refs. 10 and 11, 








R'l = 27.9 (11.0) 



Figure 1. - Contacting disk geometry (ref. 10). (Ail linear dimensions 
in cm (in. ). ) 


22 


Dimensionless measured film thickness, 


jy ' /d 



Maximum Hertz stress. 
N/m^ (psi) 

— o 1.04x109(150 000) 

■— D 1.38x109(200 000) 

— c) 1.72x109 (250 000) 

...zw 2.07x109 (300 000) 

— o 2.41x109 (350 000) 


Maximum Hertz stress, 


-o 1. 04x10' (150 000) 

_□ 1.38x109(200 000) 

,_o 1.72x109(250 000 ) 

2.07x109(300 000 ) 

..o 2.24x109(325 000 ) 





2 4 6 8 10 

(a) Type II ester lubricant. 


(b) Fluorocarbon lubricant. 


Maximum Hertz stress, 


-o 0.69x109(100 000) 

-□ 1.04x1o9 (150 000) 

-D 1. 38x109 (200 000) 

1.72x109 ( 250 000) 

JO 2.07x109 (300 000) 



Maximum Hertz stress, 
N/m^ (psi) 

- 1.04x109 (150 000) 

- 1.38x109 (200 000) 

- 1.72x109 (250 000) 

- 2.07x109(300 000) 

- 2.41x109(350 000) 


(c) Polyphenyl ether lubricant. 


Dimensionless speed-viscosity parameter, (pgU/E'R') 

(d) Synthetic paraffinic lubricant (data from ref.' 22). 


Figure 2. - Effect of maximum Hertz stress on sensitivity of minimum film thickness to changes in speed and viscosity for X-ray test 
lubricants (refs. 10 and 11). 





(a) Type Hester. 


J I I 

4 5 6 

(b) Fluorocarbon. 


7x10 


-3 



3 4 5 6 7 8x10 ' 

Dimensionless stress parameter, 

(c) Polyphenyl ether. (d) Synthetic paraffinic. 

Figure 3. - Effect of speed and viscosity on dependency of film thickness to maximum Hertz stress. 





Figure 4. - Effect of maximum Hertz stress 
on measured minimum film thickness. 



Figure 5. - Variation of EHD high-contact- 
stress factor with operating contact stress. 
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Dimensionless film thickness parameter, 


200x10'® 




temperature, Tg=422 K(300®F). 


200x10'® 




Maximum Hertz stress, N/m 


4x10® 


1 


4x10’ 


Maximum Hertz stress, psi 


(c) Fluorocarbon; disk (d) Fluorocarbon; disk 

temperature, Tg = 422 K temperature, Tg = 533 K 

(300° F). (500° F). 


Figure 6. - Comparison between predicted minimum film thickness and 
X-ray test data (refs. 10 and 11). 
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Dimensionless film thickness parameter, 








Dimensionless fiim thickness parameter, 


to 

00 


Isothermal theory, Cheny (ref. 23) 

Thermally corrected theory, Cheng (ref. 18) 

X-ray test data (from fig. 6) 



<o 


Shear heating 
KXKIxlO"^ reduction factor. 




Dimensionless stress parameter, = p^^^E' 

(a) Disk speed, 20 000 rpm. (b) Disk speed, 5000 rpm. 

Figure 7. - Comparison of X-ray measured film thickness with 
calculated isothermal and thermally corrected film thickness 
for a synthetic paraffinic oil at disk temperature Tg = 339 K 
(150° F). 


60x10'^ 



I I I I I I I I I I I 

100 200 300 400 500 600 

Disk temperature, °F 

Figure 8. - Relative fiim forming ability of test lubricants at elevated temperatures. 
Maximum Hertz stress, 2.07x10: newtons per square meter (300 000 psi); dimen- 
sionless speed factor, G = (u/Nr*®^ where N = 9.4 meters per second (370 in. /sec) 
for valid limits of 1< G < 2. 36. 
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